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ABSTRACT 

We discuss the nature of faint 6.7 /im galaxies detected with the mid-infrared camera ISOCAM on 
board the Infrared Space Observatory (ISO). The 23 hour integration on the Hawaii Deep Field SSA13 
has provided a sample of 65 sources down to 6/iJy at 6.7 /mi. For 57 sources, optical or near-infrared 
counterparts were found with a statistical method. All four Chandra sources, three SCUBA sources, 
and one VLA/FIRST source in this field were detected at 6.7 /xm with high significance. Using their 
optical to mid-infrared colors, we divided the 6.7 /jm sample into three categories: low redshift galaxies 
with past histories of rapid star formation, high redshift ancestors of these, and other star forming 
galaxies. Rapidly star forming systems at high redshifts dominate the faintest end. Spectroscopically 
calibrated photometric redshifts were derived from fits to a limited set of template SEDs. They show 
a high redshift tail in their distribution with faint (< 30/iJy) galaxies at z > 1. The 6.7 fim galaxies 
tend to have brighter K magnitudes and redder I — K colors than the blue dwarf population at 
intermediate redshifts. Stellar masses of the 6.7 /urn galaxies were estimated from their rest- frame 
near-infrared luminosities. Massive galaxies (M sta r ~ 10 11 M Q ) were found in the redshift range of 
z = 0.2-3. Epoch dependent stellar mass functions indicate a decline of massive galaxies' comoving 
space densities with redshift. Even with such a decrease, the contributions of the 6.7/xm galaxies to 
the stellar mass density in the universe are found to be comparable to those expected from UV bright 
galaxies detected in deep optical surveys. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: luminosity function, 
mass function — galaxies: stellar content — infrared: galaxies — surveys 



1. INTRODUCTION 

The evolution of galaxies can be traced in their spec- 
tral energy distributions (SEDs) as a result of the ag- 
ing and accumulation of stellar populations. Because 
more massive stars evolve faster, ultraviolet (UV) emis- 
sion becomes relatively weaker as time elapses, provided 
that no further star formation occurs. This sensitivity 
to on-going star formation can be used to trace the star 
formation history in galaxies; however, the UV light suf- 
fers from dust extinction. Thus, optical observations of 
the redshifted UV emission from distant galaxies need 
undesirable corrections if we are to deduce the evolution 
of their star forming activity. 

However, if we observe galaxies at a longer wavelength, 
this undesirable situation is improved. At rest near- 
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infrared wavelengths, much of the emission originates 
from low mass stars. Their lifetime is comparable to 
the age of the universe; thus, the effect of aging is much 
milder than that in the UV. The complicating effects of 
dust extinction are almost negligible at this wavelength. 
Most of the changes in the near- infrared SEDs are caused 
by the accumulation of stellar populations in galaxies. 
These facts assure good accuracy in estimating stellar 
masses of galaxies from their near-infrared luminosities. 

The stellar mass is one of the most fundamental 
quantities of galaxies. This motivates the execution of 
galaxy surveys in the near-infrared. For distant galax- 
ies, very sensitive surveys should be performed at a 
longer wavelength in order to detect their rest-frame 
near- infrared light. T he mid-infrared camera ISOCAM 
ijCesarskv et al. J[ 1996D on board th e Infrared Space Ob- 
servatorv (ISO. lKessler et all!996|) has achieved such re- 
quirements for the first time. 

Most of the deep ISOCAM surveys have utilized 
two broad band filters: LW2 (5 -8.5 um) at 6.7 urn 
and LW3 (12-18 urn) at 15 itm (lAltieri et all 119991: 
Flores et all ll999alU: JO liver et all 120021: iSerieant et alJ 
1997t iTaniguchi et al.1 IIQQtTT "in particular, those at 
15 urn attracted much interest l)Franceschini et aLll2001I 
120031 iGenzel fc Cesarskvll2000ft . mainly due to the dis- 
covery of a strongly evolv ing population of galaxies below 
lmJy l|Elbaz et al]ll999D . The excess in number could 
be explained by a large number of star forming galaxies 
at z ~ 1. 

The 6.7 /mi cosmological observations with ISOCAM 
resulted in generally fewer detections than at 15 /jm and 
thus attracted less attention. This is because such ob- 
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servations should sample a wavelength valley in galaxy 
SEDs between stellar and hot dust emission. The pass- 
band of the LW2 filter matches the location of the 
unidentified infrared band (UIB) emission in the local 
universe. Such dus t emission shifts aw ay from the pass- 
band with redshift ijAussel et al.ll999tL At high redshifts 
z > 1, the major contribution to the LW2 band be- 
comes stellar emission in galaxies. However, such emis- 
sion from distant galaxies requires extremely high sen- 
sitivities. Based on the importance of detecting near- 
infrared stellar emission from high redshift galaxies for 
estimating their stellar masses, we conducted a very 
deep 6.7 urn surv ey in the Hawaii Deep Field SSA13 
ijSato et al.ll2003|) . One sigma sensitivity reached 3//Jy, 
which is smaller than the p ublished values of 7 /iJy in 
iTaniguchi et all 1119971) and inlAussel et alJ l|1999|) . Most 
recently HM^tcalfe et al.l l)2003|) reach a depth similar to 
ours using a massive cluster lens. In this paper, we dis- 
cuss insights deduced from the faint 6.7 fim galaxies de- 
tected in the SSA13 field. 

After describing the 6.7 /jm sample in Sect. [21 we dis- 
cuss the identification of multi-wavelength counterparts 
(Sect. |3J). In Sect. 21 the nature of the identified galax- 
ies is examined w ith the GRASIL star/dust SED model 
ijSilva et al.lll998() . Their stellar masses are obtained us- 
ing the rest-frame near-infrared luminosities, and then we 
discuss the evolution of the stellar mass function and the 
stellar mass density in the universe (Sect. 03). Finally, we 
present a discussion (Sect. and conclusions (Sect. d>. 
Throughout this paper, we assume a flat universe with 
Q M = 0.3, Q A = 0.7, and H = 65 kms" 1 Mpc -1 . All 
optical and near-infrared magnitudes in this paper are in 
the Vega system. 

2. THE 6.7 Atm SAMPLE 

A deep mid-infrared survey has been conducted with 
the ISOCAM array on board the Infrared Space Obser- 
vatory (ISO). The broadband filter LW2 (5-8.5 fim) with 
the reference wavelength of 6.7 /im was used to image a 
high galactic latitude region in the Hawaii Deep Field 
SSA13. Many raster observations totaling an observ- 
ing time of 23 hours were combined to give a map with 
a nominal areal coverage of 16 arcmin 2 with the beam 
FWHM of 7.2 arcsec. D etails of the observations are 
given in lSato et all l|2003D . 

With a non-uniform noise distribution over the map, 
source detections were performed using a signal-to-noise 
ratio (S/N) map. The noise map was created from the 
standard deviations of the co-added raster images at each 
pixel. Detection parameters were determined by com- 
paring numbers of detected sources at the positive and 
negative sides of the map (N pos and N neg , respectively) 
with the same detection parameter set. We have cho- 
sen the parameter set giving the largest Ap OS — N ncg , 
where N pos = 65 and N neg = 12. Total 6.7 fim fluxes 
of the positive sources resulted in a range from 6 fiJy to 
170 uJv (Table EH- 

By setting a threshold that can exclude all the neg- 
ative sources, we extracted a subsample of the positive 
sources. This subsample, the primary sample, consists of 
33 sources having total fluxes larger than 12 /iJy and de- 
tection S/Ns larger than 4.3 (Table lA"T)l . Because there 
should be no effects from spurious detections in the pri- 
mary sample, our main results will be deduced based on 



the primary sample. No corrections for the contamina- 
tion of fake sources are necessary for the primary sample 
to derive integrated quantities such as stellar mass func- 
tions and stellar mass densities discu ssed in Sect . 03 In 
fact, galaxy number counts shown in ISato et alJ (2003) 
were obtained only with the primary sample. 

The remaining 32 positive sources are put into the sup- 
plementary sample. For the number of negative sources 
(-^neg = 12) having comparably low significance values, 
the number of true sources in the supplementary sam- 
ple could be 20. However, we expect more sources in 
the supplementary sample are real. This is because the 
adopted image processing could produce negative ghosts 
around bright sources, and at le ast three of the n egative 
sources were identified as such (Sat o et aLll2003|) . Gen- 
uine sources in the supplementary sample should provide 
meaningful information on the 6.7 /im sources, especially 
at the faintest flux levels. Some of them will be men- 
tioned in the following sections with a caveat that there 
will be some effects from spurious sources. 

In Table [ATI sources in both samples are listed. Here 
we also show the 12 negative sources as the negative sam- 
ple. All are to be examined with the source identification 
procedure in the next section. 

3. SOURCE IDENTIFICATION 

Counterparts of the 6.7 fim sources were searched at 
multiple wavelengths, from the X-ray to the radio. We 
used two identification methods, the probability ratio 
method in the optical and near-infrared, and the near- 
est neighbor search at X-ray, submillimeter, and radio 
wavele ngth s. The results are summarized in two tables, 
Tables IA1I and IA2I respectively. 

3.1. Optical and Near- Infrared Identifications 
3.1.1. The identification procedure 

Optical and near-infrared counterparts of the 6.7 fim 
sources were sought in the AT, /, and B bands. 
These three band data were taken from the ground 
l|Cowie et alJT9 9fi). but we utiliz ed a g5T/WFPC2 cat- 
alogu e in the isi4 band as well ijCowie. Hu. fc Songailal 
1995). The catalogue limits were estimated as K = 19.9, 
J 8 i4 = 24.7, / = 23.7, and B = 25.5 (3 a). Photometric 
uncertainties were derived by summing relative and ab- 
solute uncertainties quadratically. Absolute photometric 
uncertainties were estimated to be 0.1 magnitud e for cor- 
rected aperture magnitudes ijCowie et al.lll994) . We do 
not distinguish between igi4 and I magnitudes in the fol- 
lowing, because these magnitudes for the same sources 
in our sample do not show significant differences. At 
these magnitude limits, surface densities become higher 
than that of the 6.7 /im sources. We then introduced the 
following two probabilities to evaluate whether a candi- 
date counterpart of magnitude mi a t distance n should 
be considered as a true association l|Flores et al.l ll999a; 
iMann et allll997|) . 

Even for a true association between sources at two 
wavelengths, there is expected to be a certain amount 
of displacement in their positions. If we can neglect dif- 
ferences in light profiles at the two wavelengths, a major 
cause for the displacement will be measurement errors 
at both wavelengths. Based on smaller beam sizes in the 
optical and near-infrared (~ 1 arcsec or less) , we only 
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took account of measurement errors in the 6.7 /xm coor- 
dinates. With an assumption that the errors follow a 
Gaussian distribution, the probability for a true associa- 
tion P t is defined as 

P t = l-^-J exp(--)dx, (1) 

where a(= ri/r ) is the displacement of a potential coun- 
terpart normalized to the one sigma positional accuracy 
ro of the 6.7 /mi source in question. This ro value was 
determ ined with Monte-Carlo simulations in lSato et alJ 
(20.Q3). 

As long as the displacement has a finite value r\, we 
can not exclude the possibility that an irrelevant source 
will be found at a distance smaller than . Such a chance 
event should follow Poisson statistics for an area of irri 2 . 
Then the probability for a chance association P c is de- 
rived as 

P c = 1 - exp(-7rri 2 jV(< mi)), (2) 

where N(< mi) is the surface density of optical or near- 
infrared sources with magnitudes brighter than mi, the 
magnitude of a potential counterpart. Be cause only 
two s tars are expected in this 6.7 /im sample (Sato et al. 

m) was deri ved by integrat i ng di fferential 
galaxy counts shown in iMetcalfe et al.1 l)2001|) and in 
Maihara eTafl (|200lj) . 

These two probabilities, Pt and P c , were calculated for 
all potential counterparts up to a distance of 3ro from 
each 6.7/xm position. The large search radius was neces- 
sary to take ac count of the large uncertainty in the deter- 
mination of ro ( Sat o et al.l2003|) . We then selected an op- 
tical or near-infrared counterpart having the highest ratio 
of Pt/P c - If there is no candidate source with Pt/P c > 1, 
the 6.7 /im source was regarded as unidentified. Assum- 
ing that SEDs of the 6.7 /an sources are smoothly con- 
nected, we performed the identification procedure start- 
ing from the nearest wavelength. The adopted K and 
^814 band catalogues do not cover the full 6.7 fim field. 
Taking also account of the deeper limit for the Isu band 
catalogue, the identification scheme was then set from 
the K, igi4, /, to B band. If the identification at a par- 
ticular band resulted in a success, counterparts at the 
shorter wavelengths were found with a nearest neighbor 
search. A search up to 1.4 arcsec for it turned out to be 
enough. 

3.1.2. Identification results 

Of 65 positive sources, 57 were identified in the opti- 
cal or near-infrared. Fig. IA1I shows their appearances at 
the wavelength where the probability calculations were 
performed. In Table IA"T1 identification band (ID), dis- 
tance (ri), two probabilities (Pt and P c ), and their ratio 
are listed. The coordinates of the identified sources have 
been updated with those of the optical or near-infrared 
counterparts. We also show the second highest ratio of 

(2) 

Pt/Pc (Pt/Pc if other candidates were found in the 
search radius. There are relatively few cases of ambigu- 
ous identifications (i.e., where the second largest Pt/P c is 
of the same order of magnitude as the largest). In partic- 
ular, in the primary sample, only one source (#59) seems 
to fall in this category. In the following, the identifica- 
tion results for the primary and supplementary samples 
are examined separately. 



For the primary sample, we found that many sources 
were identified in the K band. All the sources that 
were identified at shorter wave lengt hs have large Pt/P c 
values. As indicated in Table IA3I they either lack K 
data ( #30, #39, and #6 2) or have blue colors (#44 and 
#64). ISato et alJ lj2002|) identified the source #40 with 
a submillimeter source with a hard X-ray counterpart. 
This identified source has a smaller significance value 
(Pt/Pc ~ 0.6) than the nominal threshold (P t /P c = 1) at 
the K band. We regard it as likely that the asy mmetric 
profile of this 6.7 /xm source l)Sato et al.ll2002|) affected 
the probability calculation. The only unidentified so urce 
#27, having a high S/N=10, lacks K data fTable lA3)l 
and has a very red mid-infrared to optical color. 

The lower significance of the supplementary sample 
leads to poorer positional accuracy (ro ~ 2-3 arcsec), 
which made the identification more difficult. Even in 
such circumstances, 78% (25/32) of the supplementary 
sample were identified. Because the supplementary and 
negative samples have the same S/N level, the identi- 
fication significance for the supplementary sample can 
be checked by a comparison with the negative sample. 
No source in the negative sample was identified in the 
K band. Thus the identifications of the supplementary 
sample in the K band are expected to be robust. Three 
sources in the negative sample were identified in the / or 
isi4 band. However, most of the / or igi4 band identifi- 
cations in the supplementary sample have better statis- 
tics, i.e, larger Pt/P c values than those for the I or Jgi4 
band identifications in the negative sample. Here / and 
/8i4 counterparts should be treated separately because 
of the difference in the catalogue depth. The deeper /si4 
sources generally have smaller Pt/P c values due to the 
larger surface density of sources. The 7si4 counterparts 
for sources #21 and #24 have larger Pt/P c values than 
that of the isi4 counterpart for the negative source N7. 
The I counterparts for sources #12, #29, #34, #46, 
#58, and #60 have larger Pt/P c values than that of the 
I counterpart for the negative source N9, which has a 
Pt/Pc value larger than N10. The remaining / or Ism 
band identifications are for sources #2, #31, #43, and 
#52, though most of them have slightly smaller Pt/P c 
values than that of their corresponding negative source: 
N7 or N9. Thus, we expect the number of erroneous 
identifications in the supplementary sample will be quite 
small, four at most. The number of unidentified sources 
in the supplementary sample (7) is comparable to that 
in the negative sample (9). This indicates that some 
of the unidentified sources in the supplementary sample 
could be spurious, especially those detected in noisy re- 
gions of the map (e.g. sources #7 and #37). However, 
some could be very red sources like the source #27 in the 
p rimary sample. 

iFlores et all (|1999a|) give 71 and P c values for the iden- 
tification of their 6.7 /xm sources. Their values are gen- 
erally larger than ours; the median r\ is 4.2 arcsec and 
P c exceeds 0.3 in some cases. We think that this may 
be partly due to their neglect of distortion corrections 
in data processing, which could be as large as one 6 arc- 
sec pixel in size. They adopted 1.5 arcsec map pixels and 
used only 3-5 sources to determine shifts between images. 
The median FWHM of 6.7 /im sources on their final map 
was 11 arcsec, while ours was 7.2 arcsec. We took ac- 
count of distortion corrections and used 0.6 arcsec map 
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pixels and 19 referen ce sources for the image registra- 
tion ijSato et al]l2003th The difference mentioned above 
might not explain all of the difference in their identifi- 
cation results. Rather, even the sources id entified with 
the least significance in lFlores et alJ l)1999aj) have 6.7 /jm 
fluxes larger than ours and correspondingly optical mag- 
nitudes brighter than ours. 

3.2. X-ray, Submillimeter, and Radio Identifications 

At the flux levels of the adopted X-ray, submillimeter, 
and radio catalogues, the surface densities are lower than 
that of our mid-infrared sample. Thus, source identifica- 
tion at these wavelengths can be achieved using a simple 
method, a nearest neighbor search. 

With the deep Chandra observations of SSA13 
(Mus hotzkv eTalT l200Ch . four X -ray sources were de- 
tected within the area of the ISOCAM survey. Three of 
them were detected in the hard X-ray band (2-10 keV) 
with fluxes larger than 4 x 10~ 15 ergs -1 cm -2 , while two 
of them were detected in the soft X-ray band (0.5-2 keV) 
with fluxes larger than 8 x 10~ 16 ergs _1 cm~ 2 ( Table |A"2|) . 
Within l.larcsec of these four X-ray sources, compara- 
ble to the positional accuracy of the Chandra satellite, 
we found 6.7 /im counterparts #14, #17, #40 and #62. 
Here we have used the 6.7/xm coordinates updated with 
the optical and near-infrared identifications above. 

The Submillimctrc Common-User Bolometer Array 
(SCUBA) on the 15 m James Clerk Maxwell Telescope 
(JCMT) detected three 850 /xm sources d own to 2mJy 
in th e deep SCUBA survey of SSA13 lIBarerer et all 
1999) . Taking account of the very broad beam at 850 /im 
(15arcsec FWHM), all the three submillimeter sources 
were identified with 6.7 /xm sour ces #28, #40, an d #57 
at distances of up to 7.7arcsec llSato et al .112001 . Note 
that one of the submillimeter sources (#40) was also de- 
tected in the hard X-ray. 

The VL A/FIRST survey at 1.4 GHz detected one 
3 mJy source in the 6.7/xm SSA13 map l)White et al.1 
Il997|i . This source is also l isted in the VLA/NVSS cat- 
alog ue l|Condon et a.l .1 119981. The 6.7 /im counterpart is 
assigned to the source #20 at a distance of 0.2 arcsec. 

All these X-ray, submillimeter, and radio sources are 
in the primary sample. Except for one optically bright 
hard X-ray source #62, the others have 6.7/im fluxes in 
a range of 10-30 /xJy. 

4. THE NATURE OF FAINT 6.7/im GALAXIES 

We have identified two stars in the 6.7 Mm sample us - 
ing their image profiles in the optical ijSato et alJ l2003'). 
Excluding these stars from the primary sample, we now 
discuss the nature of the remaining sources, all of which 
are assumed to be galaxies. We utilize the G RASIL 
star/dust galaxy SED model (|Silva et al.lll998j) . which 
is explained in some detail in Appendix lAl 

4.1. Colors 

With the K, I, and B band photometry, we examined 
the distributions of faint 6.7 um galaxies in various color- 
color plots. In one plot (Fig. EH), we identified a distinct 
population of faint 6.7 /im galaxies having red B — I col- 
ors and low ratios of /j, (6.7/im) //„ (if). To derive ratios 
with 6.7/im fluxes, we used zero-point fluxes of 645 Jy, 
2408 Jy, and 3974 Jy for K, I, and B magnitudes, respec- 
tively. Solid circles with error bars show sources in the 



primary sample. Empty circles are for the supplemen- 
tary sources. Color limits using 3 a magnitude limits or 
real measurements for saturated sources or for sources at 
the image boundaries are indicated with arrows. Some 
sources are not shown due to lack of information, such 
as no data in the K band, or no detections in both I and 
B bands. The concentration of galaxies at the upper-left 
part of the panel is distinct from the rest of the sample by 
more than the relatively large errors in f v {<o.7^,m)/f v {K) 
ratios. We assigned an identification flag of type I for 
galaxies in this concentration (Table IA3|) . The type I 
galaxies are marked with squares and can be separated 
from the rest of the galaxies with the dot-dashed line in 
the plot. 

On the same panel, we overlaid some model predictions 
using the GRASIL SEDs (Appendix EJ). The evolving 
SEDs for three z = Hubble types - E, Sa, and Sc galax- 
ies - are shown with dashed, dotted, and solid lines, re- 
spectively. For each Hubble type, we assumed three for- 
mation redshifts of Zf — 2, 3, and 10 (thick, medium, and 
thin lines, respectively). The locations of z = galaxies 
are shown with solid triangles, while z — 1,2, and 3 with 
empty triangles (for z < Zf). We find that type I galax- 
ies follow lines for the evolving E galaxies at intermediate 
redshifts of z — 0-1, regardless of their assumed forma- 
tion redshifts. The evolving Sa galaxy predictions also 
share the same region, but higher formation redshifts are 
preferred. In our adopted cosmology, the ages of galax- 
ies with Zf = 2, 3, and 10 at z = 1 become 2.5 Gyr, 
3.6 Gyr, and 5.3 Gyr, respectively. This indicates that 
type I galaxies would be old and matured systems whose 
stellar contents are already in place as a consequ ence of 
vigorous star formation a long time ago (see Fig. lAl4| . 

The other galaxies in this plot follow the GRASIL pre- 
dictions; however, almost nothing can be drawn from 
this figure because of the degeneracy of the model loci. 
Both rapidly star-forming galaxies at high redshifts and 
quiescently star-forming galaxies show high ratios of 
/„ (6.7/im)//^ (K) and blue B — I colors. The color de- 
generacy in this part of the plot can be disentangled to 
some degree by changing color combinations. In Fig. lA3l 
we show / !/ (6.7/im)// 1/ (7) ratios and B — K colors. In 
this plot, the evolving E or Sa galaxies can be red in 
both colors, while the evolving Sc galaxies occupy the 
lower left part of the panel with blue B — K colors and 
low ratios of /„(6.7/xm)//„(J). The evolving E and Sa 
galaxies with lower formation redshifts (zf < 3) can have 
such blue colors at their forming stage (z ~ Zf). Thus, 
these blue colors can be an indicator of on-going star 
formation. 

The distribution of the 6.7/xm galaxies in this plot is 
rather smooth, except for the type I galaxies marked with 
squares. However, we here introduced a separation line 
(dot-dashed line) to extract many sources not detected 
in the B band. These can be explained as high redshift 
galaxies (z > 1) with rapid star forming activities at the 
past (evolving E or Sa galaxies). They represent post- 
starburst galaxies at high redshift. Their properties are 
consistent with those of ancestors of type I galaxies. We 
here categorize them as type II and mark them with dia- 
monds. The remaining galaxies are assigned as type III. 
They are on-going star formers at z = 0-2, a combina- 
tion of mildly star forming galaxies and vigorously star 
forming galaxies with low formation redshifts. 
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The division of the sources into types I, II, and III 
is made in color-color plots with photometry in all four 
bands - 6.7/j,m, K, I and B. We find that this di- 
vision can also be applied to the color-color plot of 
f v (6.7fj,m)/f v (I) vs f u (6.7^im)/f u (B). In this plot, we 
can assign type identifiers to all the 6.7 /jm sources, 
with the help of the distribution of the sources whose 
types have been defined in the previous two plots. Even 
unidentified sources could have a nominal type of II by 
their red colors due to their flux limits at the I and B 
bands. For the primary sample, the number ratio for 
the 31 galaxies becomes type I : II : III = 11 : 10 : 10 
inc ludin g the unidentified source #27 with type II (Ta- 
ble IA3|I . This almost even distribution skews toward the 
dominance of type II in the supplementary sample. For 
the 32 sources in this sample, the ratio becomes type 
I : II : III = 5 : 23 : 4. The dominance of type II in this 
fainter sample is still valid even if the seven unidentified 
sources assigned type II are excluded. 

4.2. Redshifts 

In order to obtain absolute quantities such as stel- 
lar masses for the 6.7 /im galaxies, we need their red- 
shifts. For the 55 identified 6.7/im galaxies, we do have 
21 spectroscopic redshifts z S n, of which 18 are for the 
primary sample l|Barger et alJ 120011 iCowie et alJ Il996t 
iSongaila et alJll994T) . For others, we derived photomet- 
ric redshifts z p h in the following way. 

4.2.1. Photometric redshifts 

Photometric redshifts were estimated by minimizing 
X 2 values between photometric measurements and model 
estimates. This is basicall y the same method adopted 
in the the hyper z code bv IBolzonella. Miralles. fc Pellol 
( 2000). We used 6.7 /i m, K , I, and B band fluxes or t heir 
limits listed in Table IA3I For sources in Table IA2I we 
tried to utilize their flux values or their limits at the X- 
ray, submillimcter, and radio wavelengths as well. Model 
estimates were obtained using two sets of the star/dust 
SEDs in the GRASIL library; the evolving SEDs and 
the local ones (Appendix 0. We also added two SEDs 
expanded to the X-ray wavelengths; a NGC 6240 SED 
compiled by lHasingerl p000|) and a GRASIL Arp220 
SED complemented with X-ray an d radio observations 
ICarilli fc YurJ f2000: I was aw 11999T) . 

The x 2 minimization was executed SED by SED in our 
code. The SEDs used have a long wavelength baseline; 
the GRASIL local SEDs from UV to submillimeter, the 
GRASIL evolving SEDs from UV to radio, and the two 
additional SEDs from X-ray to radio. They are likely 
to have significance variations at different wavelengths; 
however, we neglected such errors in the SEDs. We thus 
considered noise only in the absolute photometric mea- 
surements. For each SED, we calculated x 2 values at 100 
redshifts from z = 0.1 to 10. The interval of the redshift 
grid was constant in log(l + z). Once the redshift giv- 
ing the minimum x 2 value was derived, we repeated the 
calculation around this redshift with a ten times finer 
redshift resolution. After completing the x 2 minimiza- 
tion for all the SEDs, we determined a photometric red- 
shift giving the minimum x 2 value among the SEDs. We 
also found that in most cases, the redshift giving the 
minimum x 2 value for each SED is consistent with the 
photometric redshift determined by the global minimum 



in x 2 ■ It should be noted that the final photometric red- 
shifts were all determined with the UV-to-submillimeter 
SEDs because the SEDs expanded to the radio and/or 
X-ray resulted in high x 2 values. For the comparison of 
X 2 values among SEDs with different wavelength spans, 
we took proper account of both the values of x 2 and the 
number of the detected data points used to calculate x 2 ■ 
In these calculations, we did not require that the age of 
the SED template be less than the age of the universe 
at a given redshift. Previous works on photometric red- 
shifts have given good results with loca l templates at 
high redshifts (see e.g.. lHogg et alJ ll998). We therefore 
tried to fit local SEDs even at z = 10. Old evolving SEDs 
were also treated as local SEDs. This assures a possible 
range of SED variations. Actually, we found that the 
results of this approach admitting the age-inconsistency 
in the fits gave the best results in the comparison with 
spectroscopic redshifts explained below. 

We sh ow some fitting results at different redshifts in 
Fig. lA4l Ph otome tric redshifts and th eir 9 % confidence 
limits l|Avnilll97ql are listed in Table lATil The names of 
the best-fit SEDs, their x 2 values, and some ancillary in- 
formation of the fits are also shown. It should be noted 
that many other SEDs gave quite similar redshift esti- 
mates with somewhat larger x 2 values. For source #38, 
which has the largest x 2 value in the sample, 12 SEDs 
gave z = 0.74-0.98 with x 2 = 9.9-15. On the other hand, 
many sources have x 2 = 0.0, though they are truncated 
values. This indicates overestimates of our photometric 
errors. This is likely in our case due to quadratic sum- 
mation of error components, each of which was difficult 
to determine independently Our data points are very 
sparse in wavelength. In a sense, we might have used 
too many sets of SEDs for such data sets. It can be 
said that we have allowed the SEDs to shift freely in red- 
shift space to take account of the symmetric distribution 
of photometric errors. We have succeeded in eliminat- 
ing outlier fits; however, some fits could have fallen at 
rather artificial x 2 minima. Thus, the confidential inter- 
vals are more reliable than photometric redshift values 
themselves. It should be noted that the photometric red- 
shift determination technique basically utilizes the char- 
acteristic spectral features common to all SEDs. There 
is some information in the fit; however, the very particu- 
lar SED names identified as giving the best photometric 
redshifts do not exclude other SED types. 

The reliability of these redshift estimates can be ad- 
dressed by computing photometric redshifts for sources 
with spectroscopic redshifts. The results are shown in 
Fig. IA5I Of 21 sources with spectroscopic redshifts, 
10 photometric redshifts were derived with the age- 
inconsistent evolving SEDs, 6 with local templates, and 
5 with age-consistent SEDs. This result assures the 
usefulness of the admittance of age-inconsistency in the 
fits. Actually, most of the photometric estimates follow 
the identity relation with their spectroscopic measure- 
ments. The nominal dispersion was 0.2 dex in logz (dot- 
ted lines). Only four sources, #62, #15, #29, and #16 
with spectroscopic redshifts of z ~ 0. 2, 0.3, 0.6, and 1.2, 
respectively, gave photometric redshifts larger by more 
than one sigma. However, their large 90% confidence 
intervals indicate that they are essentially problematic 
data sets. For example, source #62 lacks any flux con- 
straint in the K band. 
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The four photometric bands, B, I, K, and 6.7 /i m, ar e 
widely and evenly distributed in wavelength (Fig. IA4|) . 
The determination of photometric redshifts with these 
data sets appears to depend on a characteristic spec- 
tral peak around 1.6 /im corresponding to the H~ opac- 
ity minimum for the stellar continuum emission. The use 
of this 1.6 /xm bump for photometric redshifts with pho- 
to metric data set s without wavelength gaps is described 
in iSawickll (|2002^ . In wider wavelengths, the spectral 
slope is generally monotonic below 1.6 /im, while there 
could be UIB emission longward of the 1.6 /im peak until 
6.7 /zm (i.e. 3.3 /jm and 6.2 /im). There is a possibil- 
ity that some UIB emission was erroneously treated as a 
shorter wavelength feature. Misidentification of 6.2 /im 
emission with a 3.3 /tm or 1.6 /im feature gives a wrong 
redshift of z ~ 1 or 3 and that of 3.3 /tm emission with 
a 1.6/tm feature gives z ~ 1. Unfortunately, the 6.7/tm 
photometric band is very broad (5-8.5 /tm) and the gaps 
between the four photometric bands are somewhat large. 
Some sources are even lacking some of the four flux val- 
ues. This all c ould give a reason for the four high redshift 
outliers in Fig. lA5l The effects of z p h > z sp estimates will 
be addressed in the following sections. 

4.2.2. Fluxes and redshifts 

With redshifts determined spectroscopically or photo- 
metrically, 6.7 /im fluxes for all 55 ident ified galaxies are 
plotted as a function of redshift in Fig. IA6I Flux errors 
are indicated with vertical bars, while horizontal bars are 
90 % confidence limits for photometric redshifts. Sources 
with spectroscopic redshifts are marked with double cir- 
cles. Galaxies in the primary sample are highlighted with 
solid symbols. 

For the primary sample, fainter galaxies tend to have 
higher redshifts. This usual trend seems to weaken be- 
low a flux level of 30 /iJy, where high redshift galaxies 
start to appear regardless of their flux values. This 
first finding is a direct result of the substantial depth 
of our imaging. Some photometric redshifts could be er- 
roneously high (Sect. I4.2.1fl . but the high redshift tail 
can be recognized even within the spectroscopic primary 
sample alone. There are some z sp > 1 galaxies with 
fluxes of ~ 20/iJy. Most of the supplementary sample 
have z p h > 1. Because they are generally fainter than 
30 /iJy, this behavior itself strengthen the high redshift 
tail seen in the primary sample. Two bright supplemen- 
tary sources at z ~ 2 (#12 and #34) are detected in rel- 
atively high noise regions and lack K data. Thus, their 
significance is taken to be low. 

In this plot, model predictions are overlaid using the 
GRASIL library (Appendix EJ . Dashed, dotted, and 
solid lines are evolving E, Sa, and Sc galaxies and thick, 
medium, and thin lines are cases with formation redshifts 
of Z( — 2, 3, and 10. The distribution of the 6.7/im 
galaxies is almost bracketed by these model predictions. 
At low redshifts (z ~ 0.3), most of the 6.7/im galaxies 
are consistent with the evolving Sa and Sc galaxy mod- 
els. The steep flux-redshift slopes for the evolving Sc 
galaxies are due to the effects of strong dust emission in 
the Sc galaxies to the 6.7 /tm observing band. At high 
redshifts (z ~ 1), the distributions of both primary and 
supplementary galaxies are centered at model predictions 
for the evolving Sa galaxies. The flat flux-redshift dis- 
tribution of the 6.7 /tm galaxies at high redshifts can be 



reproduced well by a factor of a few less luminous version 
of the evolving E galaxy models. Although the GRASIL 
models cannot be arbitrarily rescaled, especially at the 
UV and far-infrared wavelengths where effects of dust 
are significant ( Appendix lA"| . these model loci at z > 1 
are primarily due to stellar emission in the rest-frame 
near-infrared, which accepts rescaling. 

We mark the locations of the X-ray, submillimeter, and 
radio sources with letters, 'X', 'S', and 'R', respectively. 
Except for one optically bright hard X-ray source (#62), 
all are within the high redshift tail. 

4.3. Comparison with a K < 20 Sample 

Some of the characteristics of the faint 6.7 /tm galaxies 
can be assessed in a comparison with a K band magni- 
tude limited sample. We here utilize the K < 20 spectro- 
scopic sample in the SSA13 field presented in lGowie et alJ 
( 1996) . The field coverage of this K sample is a few times 
larger than that of the ISOCAM map. 

First, we show a K — z dia gram both for the 6.7 /inl- 
and Jf-sel ecte d samples in Fig. lA7l Symbols are the same 
as in Fig. IA6I except K < 20 galaxies are overlaid with 
sma ll sq uares. Note that some 6.7/im galaxies plotted in 
Fig. EH are not shown here because they lack K photom- 
etry. 6.7 /tm galaxies with fainter K magnitudes tend to 
have larger redshifts. This feature is well described with 
the GRASIL evolving Sa galaxy model, and the distri- 
bution of the 6.7/im galaxies in the plot lies between the 
evolving E and Sc galaxy predictions. No significant dif- 
ference is seen between the primary and supplementary 
samples. 

In contrast, the distribution of the K band selected 
galaxies extends to a region fainter than the evolving Sc 
galaxy model predictions and contrasts with the 6.7 /tm 
sample. The 6.7 /tm sample preferentially selects a high 
redshift population at each K magnitude. At intermedi- 
ate redshifts (z = 0-1), K band emission is a good indica- 
tor of stellar masses (Sect. l5~T|) . The evolving Sc galaxy 
models at z ~ 0.4 have stellar masse s of 0.1-0. 2 M„* , 
where we adopt M* tar = 1.7 x 10 11 M ijCole et al.l2001|i . 
Thus, the z < 1 faint K galaxies should have very small 
stellar masses. 

Next , we show I — K colors as a function of redshift 
in Fig. IA8I The 6.7 /tm galaxies have red I — K colors 
especially at z > 1. This reddening trend is expected 
with the GRASIL models. However, some 6.7/im galax- 
ies exceed the red envelope of the GRASIL predictions. 
These excess red I — K colors are seen at intermediate 
redshifts for sources in the primary sample, and seem 
to be consistent with the presence of some high redshift 
(z > 1) galaxies in the supplementary sample. The very 
red I — K colors at z > 1 indicate rapid star forma- 
tion at very early epochs in the context of the GRASIL 
model. The reddest I — K colors in the 6.7/im sample 
are comparable to those of an evolving E galaxy with a 
formation redshift of Zf = 10. For a brighter 6.7 /tm sam- 
ple, the existence of red I — K sources is also reported in 
IFlores et al.lfl999^. 

Most of the 6.7/tm galaxies have I — K colors red- 
der than the evolving Sc galaxy predictions. However, 
some K band selected galaxies have blue colors such as 
I — K < 2. Their blue colors imply that young stel- 
lar populations are dominant. As indicated with the 
GRASIL predictions at z ~ Z{ , very blue I — K colors can 
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only be seen in the forming stages of galaxies. Thus, the 
blue if -selected galaxies should experience star forming 
activities less burst-like than the evolving Sc galaxies, or 
they should start to form stars in a wide range of redshift 
at z = 0.3-1. 

We find that the K band galaxies bluer than the evolv- 
ing Sc galaxy predictions have K magnitudes fainter than 
the Sc galaxy models (Fig. IA7I) . These K selected galax- 
ies having blue I — K colors and faint K magnitudes at 
intermediate redshifts (z < 1) should have small stellar 
masses at their forming stages. Faint K band selected 
samples will be contaminated with this population of 
young dwarf galaxies at low redshifts. However, faint 
6.7 /im selected samples will be free from it. It means 
that an efficient search for high redshift galaxies will 
be achieved with deep surveys at the mid-infrared, e.g., 
6.7 /im. 

5. STELLAR MASS 

5.1. Stellar Mass-to-Light Ratios 

For the conversion of observed luminosities to their 
stellar masses, we calculated stellar mass-to-light ratios 
at the observed bandpasses. These ratios are to be ap- 
plied to the observed values directly. This approach is 
different from others in literature (e.g., taking the ob- 
served K band flux for a z = 1 galaxy, converting it to 
the rest-frame K band with an assumed SED, and ap- 
plying a stellar mass-to-light ratio at the rest-frame K 
band). We eliminated such conversion to a fixed band- 
pass. 

We calculated both stellar mass M st ar and an in- 
band luminosity L that should be observed directly as 
a function of redshift. The in-band luminosities are pre- 
sented in units of the bolometric luminosity of the Sun 
(1 L^ 01 = 3.85 x 10 26 [W]). We adopted the evolving SEDs 
in the GRASIL library ( Appendix IXjl and the resulting 
stellar mass-to-light ratios M star /L ar e sh own for the B, 
I, K, and ISOCAM LW2 bands (Fig.KSjl. Dashed, dot- 
ted, dot-dashed, and solid lines represent the evolving 
E, Sa, Sb, and Sc galaxies. Formation redshifts were 
assumed to be Z{ = 5 and 10 (thick and thin lines, re- 
spectively). 

It is well known that the rest-frame K band light is 
a good indicator of stellar mass. This fact is reflected 
in the K band panel, showing very small dispersions 
among the cases for different galaxies or different for- 
mation redshifts, especially at low redshifts. At low red- 
shifts, the dispersions among different models become 
larger at shorter wavelengths, here the I and B bands. 
The large dispersions at low redshifts in the 6.7 /jm panel 
are due to the contamination from dust emission. 

Moving to larger redshifts, the situation regarding 
the dispersion among star-forming histories will change 
among the different observing bands. Very roughly, the 
effective wavelengths of the B, I, K, and LW2 bands 
are two time larger than those of the neighboring shorter 
bands. The dispersions for the J, K, and LW2 observing 
bands at z ~ 1 are very similar to those for the B, I, 
and K bands at z ~ 0. This means that 6.7 /im light will 
become a better stellar mass indicator at high redshifts. 
On the contrary, B band fluxes cannot be used to esti- 
mate stellar masses at high redshifts. They correspond 
to rest-frame UV fluxes, which are good indicators of star 



formation rates. Emission at short wavelengths such as 
UV emission is sensitive to ages of galaxies, or formation 
redshifts. Such effects are actually seen in the B band 
panel. They can be seen in the K band panel at high 
redshifts. 

5.2. Stellar Masses of Faint 6.7 Galaxies 

Stellar masses of the faint 6.7 /im galaxies were derived 
from rest-frame near-infrared light taking into acc ount 
the behavior of the stellar mass-to- light ratios in Fig. lAQI 
Namely, the conversion of luminosities to stellar masses 
was performed using 6.7/im luminosities for sources at 
z > 1 and K band luminosities for sources at z < 1. 
For low redshift galaxies with no K photometry, we used 
/ band luminosities instead. This approach minimized 
the uncertainties in stellar mass-to-light ratios. In addi- 
tion, we determined the SED to derive stellar mass-to- 
light ratios by using all the photometry data from UV to 
submillimeter via x 2 minimization to the evolving SEDs 
in the GRASIL library (Appendix Ia")| . Because older 
SEDs usually have higher stellar mass-to-light ratios, we 
set an age constraint for the fitting template. No SEDs 
older than the age of the universe at the source's redshift 
were used for the fitting. This will avoid artificially high 
estimates of stellar mass-to-light ratios at the expense 
of increasing the nominal x 2 of the best-fit. Note that 
X 2 values for sources with spectroscopic redshifts can be 
larger than those for sources with photometric redshifts. 
This is because some fraction of inappropriateness of the 
model SEDs is absorbed in the fits in determining pho- 
tometric redshifts. However, such errors are reflected in 
the uncertainties of the photometric redshifts. 

The resulting stellar masses using the best-fit stellar 
mass-to- light ratios are listed in Table EH The names of 
the SEDs used, their x 2 values, and some ancillary infor- 
mation of the fitting are also shown. In order to derive 
errors in stellar mass-to-light ratios, we s earched SED 
fits with Ax 2 = 1 from the best-fit SED l|Avnil 1 1976ft . 
For some sources, no such SED fits were found to deter- 
mine confidence limits because of the sparse grid of SEDs. 
In such cases, we adopted stellar mass-to-light ratios of 
the SED fits with the minimum A% 2 values. Even with 
such overestimated noise, we found that the photomet- 
ric errors dominate the total uncertainties in the stellar 
masses. This is because any effects on the stellar mass- 
to-light ratio values were minimized as a result of our 
hybrid conversions using 6.7 /jm and K b and lu minosi- 
ties. The actual trend can be seen in Fig. IA10I Larger 
errors in stellar masses for sources at z > 1 originate in 
larger photometric errors in their 6.7 /im fluxes. It is also 
noted that our stellar mass estimates are not likely to be 
affected by the limitation of rescaling in the GRASIL 
SED model because the rest-frame near-infrared light is 
negligibly affected by dust (Appendix IXIl . 

Many 6.7 /im galaxies have stellar masses as massive as 
10 11 Mq. The characteris tic stellar mass of local galaxies 
M* tar = 1.7 x 10 11 Mq ijCole et all l2001|) is indicated 
with a horizontal dot-dashed line. Also shown (dashed 
line) are stellar masses of the evolving E galaxy with 
a formation redshift of Zf = 10 scaled to the minimum 
6.7 /im flux of the primary sample (12/iJy). This shows 
that M* tar galaxies could be detected out to z ~ 3, if such 
galaxies were present. Some 6.7 /im sources have stellar 
masses smaller than this GRASIL prediction. They are 
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likely to have star formation histories less burst-like than 
the evolving E galaxy or lower formation redshifts of Zf < 
10. This curve itself is useful to indicate the redshift 
dependence of the detection limit especially at the high 
redshift end where the scaling from the GRASIL model is 
almost unity. However, the scaling factor becomes very 
small, ~ 0.03 (0.1M* ai ) at z = 0.2 (Appendix © . So, 
the uncertainties of the dashed line at low redshift will 
be large. Moreover, 6.7 /im fluxes start to be affected 
by dust emission at such low redshifts. Note that stellar 
masses were derived from K fluxes at z < 1. 

The distribution of the supplementary sample is sim- 
ilar to that of the primary sample. Two of the high- 
est stellar mass sources in the supplementary sample at 
z ~ 2 (sources #12 and #34) are the dubious sources 
mentioned in Sect. 14.2.21 The lowest stellar mass for 
the supplementary sample at z ~ 0.6 (source #29) falls 
somewhat far from the distribution of other sources. This 
is derived from its / band luminosity and could be some- 
what less secure. 

5.3. Epoch Dependent Stellar Mass Functions 

Here we examine the redshift dependence of the stellar 
mass function. The epoch-dependent stellar mass func- 
tions <fr(Mst, a , ri z) were estim ated with the 1/ V max method 
(e.g.. lTakeuchi et alll20 00) as 

$(M star ,2)=£-J-. (3) 

''max 

No sources in the supplementary sample were used in 
the summation in order to avoid the unwanted effects of 
spurious detections (Sect.|2J). Taking account of errors in 
photometric redshifts, we adopted three broad redshift 
bins; z = 0.2-0.5, 0.5-1.2, and 1.2-3.0. The widths of 
these bins are almost 0.4 dex in logz. The maximum 
volume Knax for each source was derived according to 

V max = £y«S(z))^- zd z, (4) 

where il(< 5) is the effective solid angle in this sur- 
vey for a so urce with a flux of S. This function is 
presented in iSato et all l)2003f) . S(z) is the flux of a 
source when placed at a redshift z. The best-fit SED 
for the stellar mass was used to derive S(z). The comov- 

ing volume element is cal culat ed with formu lae in 

ICarroll. Press, fc Turneif l|T99^ and iHogg I (fl99^) . The 

redshift integration range is determined by the criteria 
of the primary sample. z max is the smaller of (a) the 
redshift at which the source has the minimum signal-to- 
noise ratio in the primary sample, and (b) the upper end 
of the redshift bin. z m ; n is the larger of (a) the redshift 
at which the source exceeds the maximum flux in the 
primary sample, and (b) the lower end of the redshift 
bin. 

We computed stellar mass functions both for the spec- 
troscopic sample and for the combined sample of gal axie s 
with spectroscopic and photometric redshifts ( Table lA"4|l . 
Stellar mass functions for the spectroscopic sample are 
listed to indicate the strict lower limits. One sigma er- 
rors are also shown, estimated from Poisson statistics 
and uncertainties in Vmax- The numbers in parentheses 
are the number of galaxies that were used in the calcula- 
tion. Because of the small sample size, we only used four 



broad stellar mass bins. Some values are missing due to 
the limited survey volume (at the high stellar mass end 
at low redshift) and limited sensitivity (at the low stel- 
lar mass end at high redshift). With such restrictions, 
there are few overlaps among the redshift bins for a fixed 
stellar mass bin. In each of such overlapping stellar mass 
bin, we find that comoving space densities are decreasing 
with redshift. For a stellar mass range for typical local 
galaxies (log(M sta r[/i6 5 2 M ©D = 10.95-11.35), the comov- 
ing space density at z = 1.2-3.0 becomes 10 % (—1.0 dex) 
of that at z = 0.2-0.5. 

These stellar mass funct i on estimates are compared 
with local ones in Fig. IA11I iCole et aTl l)2001l) derived a 
local stellar mass function with a large sample of matched 
2MASS-2dFGRS galaxies. Their stepwise maximum like- 
lihood estimates are shown with diamonds in the upper- 
left panel. The assumed initial mass function (IMF) was 
a Salpeter-type l)Salpeterlll955|) . The error bars are also 
shown. The minimum error is obtained almost at a char- 
acteristic stellar mass of M* tar = 1.7 X 1O 11 M0, which 
is converted to our cosmology. The redshift distribu- 
tion of the 2MASS-2dFGRS galaxies is largely confined 
to z = 0-0.2 with a peak around z — 0.05. This local 
stellar mass function is represented with dashed lines in 
higher redshift panels for the comparison with our esti- 
mates. Solid and double circles show our stellar mass 
functions derived with the combined and spectroscopic 
samples, respectively. The values for the spectroscopic 
sample are shifted leftward by 0.05 dex. Vertical error 
bars show Poisson noise. Horizontal error bars are mean 
fractional errors in stellar masses. 

In the z = 0.2-0.5 bin, our stellar mass function es- 
timates are consistent with the local ones within one 
sigma. At higher redshifts, stellar mass functions es- 
timated with 6.7 ^m galaxies start to show lower values 
than the local ones. At redshifts of z = 0.5-1.2, the devi- 
ations are somewhat marginal. The deviations appear to 
be larger at larger stellar masses. At the highest redshift 
bin z = 1.2-3.0, the deviations from the local sample are 
more than one sigma. These deviations could be much 
larger, if we took account of our tendency to overestimate 
photometric redshifts fSect. l4~2~T|l . Cosmological surface 
brightness dimming could result in a smaller number of 
galaxies at high redshifts. However, our sample has been 
selected from a 6.7 [im map with a broad beam, corre- 
sponding to a very deep surface brightness sensitivity 
(0.02/iJy arcsec -2 at la). 

The decrease in the comoving space density of massive 
stellar systems at high r edshifts has also b een seen in 
K band selected data bv lDrorv et alJ l)2001l) . They pre- 
sented integrate d stell ar mass functions above three mass 
thresholds fFig. IA12|l . which are alm ost comparable to 
our largest three mass bins (Table lA"4jl . Their integrated 
stellar mass functions (squares) for log(M s t ar [/i^5 2 M ]) > 

11.07 and log(M sta r[/i65 2 M ©D > 1L37 show declines of 
—0.5 dex and —0.8 dex from z = 0.5 to z = 1.1. Our in- 
tegrated stellar mass functions (circles) have a good con- 
sistency in the overlapping redshift range. For the low- 
est mass threshold panel, our highest redshift bin values 
are lower limits, b ecause of the limit ed sensitivity (Ta- 
ble EU). Note that lDrorv et all l)2001j) adopted the max- 
imum mass-to-light ratios as suming the age of the uni- 
verse. The local values using IQole et al.l IpOT^ 's stellar 
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mass function are marked as references (diam onds). Our 
values at z — 0.3 are somewhat larger than lCole et alJ 

(2001) 's values. This difference may result from cosmic 
variance, because our surveyed v olume at su c h low red- 
shift is much smaller than that of iCole et alJ l)200lD . 

5.4. Stellar Mass Density in the Universe 

Using the stellar mass estimates for the 6.7 /jm galax- 
ies, we derived their contributions to the stellar mass 
density in the universe. Epoch-dependent stellar mass 
densities p s t&r(z) are obtained as 

Pstar(» = TT^' (5) 

''max 

utilizing the Umax values calculated in Sect. 15.31 Here 
again, we only used the primary sample for the sum- 
mation. Our derived stellar mass densities for the three 
redshift bins are direct sums. Our stellar mass func- 
tion estimates have narrow mass ranges and large errors 
(Fig. lAllfl . They could not be used to derive charac- 
teristic masses and low-mass slopes, which are needed 
to properly integrate stellar mass functions over the full 
mass range. Thus, we decided to use direct sums in or- 
der to avoid large uncertainties in corrections to fully 
integrated values. Our direc tly summed stellar mass 
densities are shown in Table. IA4I Their errors are es- 
timated from Poisson noise and uncertainties in stellar 
mass and U max . The numbers of sources used in each 
summation are indicated in parentheses. Stellar mass 
density becomes smaller at higher redshifts. It should be 
noted that the 6.7 /im galaxies that contribute these val- 
ues come from different stellar mass ranges, as indicated 
in the left columns. 

The contributions of our 6.7 /im galaxies to the stel- 
lar mass density in the universe are shown as a func- 
tion of redshift in Fig. IA13I Solid and double cir- 
cles were estimated from the combined and spectro- 
scopic samples, respectively. Double circles are shifted 
slightly to lower redshifts. Horizontal bars show 
bin widths and vertical bars mark one sigma uncer- 
tainties. We overlay several stellar mass densitie s 
in the literatu re, as triangles (IGiallongo et alJ fl998). 
empty circles < | Brinchma nn'j£F^hT|2QnoTr a diamond 
llCc-le ei; al|2QQl l). X m arks l|CoheTl2002lUemptv squares 
(Dicki nson et alJl2003|) . and solid squares (Fontan a et al.l 
2003). 

The local value bv lCole et all IpfflOlh was obtained from 
their stellar mass function, which is deduced assuming 
a Salpeter IMF. The value is derived for the full mass 
range by i ntegrating the Schechte r fit to the stellar mass 
function. Dic kinson et alJ {2003) also adopted fully in- 
tegrated values over the full mass range. They first de- 
rived luminosity densities by integrating the Schechter 
fit to their rest-frame _B-band luminosity function at 
each redshift, and then converted them with their mean 
.B-band mass- t o-light ratio into stellar mass de nsities. 
iFontana et alJ l)2003[) followed this method. iCohenl 

(2002) used the integration of the Schechter fits to the lu- 
minosity functions, but for a restricted range from 10 L* 
to 1/20 L* . These quasi- integrated values for her if-band 
luminosity functions were converted with a fixed z = 
stellar mass-to-light ratio of M stal /LK=0-8 in solar units. 
Note that this M star /LK ratio has a different meaning 



from our M star /L ratio for the K band in units normal- 
ized to the bolomctric s olar luminosity (Sect . |Kjll . which 
is shown in Fig. |X9] iBrinchmann fc Ellis! lj2000j) also 
adopted the value of M star /ivA'=0.8 to derive incomplete- 
ness corrections. Their values are quasi-full integrations 
for a limited m ass r a nge of 10.5 < lo^jM^A h^ M Q ]) < 
11.6 (cf. Table EU). IGiallongo et all (|1998l) 's values are 
direct sums for 20 < R' < 25 galaxies around a z = 4.7 
quasar. Their values are multiplied by two to take ac- 
count of their adopted Miller-Scalo IMF. 

Some authors assumed a Schechter form; however, the 
shape of stellar mass functions is not determined yet. We 
therefore did not applied any corrections to these values 
derived for the different mass ranges. We made crude 
corrections only for the different cosmologies. Although 
the reported values are basically binned values, we took 
each of them as a single representative value at each red- 
shift with a fixed redshift range. We did not consider ef- 
fects from sources near the boundari es of the bins or with 
differe nt z m i n and z max values. For Brinchm ann fc Ellis! 
(2000), we used midpoints in logz f or their bin bo und- 
aries. We used her mean z values for ICohenl l)2002j) . 

At lower redshifts (z < 0.5), almost all the estimates 
are consistent and very near to the local value. This sug- 
gests that surveys cited here have succeeded in detecting 
most of the stellar masses in galaxies at these redshifts. 
At a higher redshift (z ~ i) ; differences between the 
results become somewhat larger. The use of a fixed stel- 
lar mass-to-ratio of M^t/ Lk = 0.8 in lCohenl l|2002fl and 
IBrinchmann fc Ellisl l)2000f > could explain so me of the val- 
ues in excess of ours at z ~ 1. According to lDrorv et alJ 
l)2001[K stellar mass-to-ratio M sta , Y /LK is a decreasing 
function of redshift; 0.99, 0.88, 0.73, and 0.65 for z = 0.5, 
0.7, 0.9, and 1.1, respectively. They assumed the maxi- 
mum ages for the passively evolving galaxies; thus, the 
decline in M stal /LK could be steeper. Note that the de- 
clines in the M stal /L ratios at the K band in Fig. IA9I 
come from the evolution of the galaxies but also from 
the effects of the shifting observing band. Direct sums 
of the detected sources only are likely to miss low mass 
galaxies below the detection limits. Therefore, our spec- 
troscopic points, which are derived from our sub-samples, 
should be regard ed as strict lower limits . At even higher 
redshifts (z > 2). IGiallongo et all l)1998l) 's values appear 
to have a milder decline, especially compared with those 
of iDickinson et aD ()2003|) . This might related to their 
use of the dif f erent IMF in the estimation. However, 
IFontana et alJ l)2003fl presented high stellar mass densi- 
ties at these redshifts. This suggests that we need to sur- 
vey much larger area to determine the true stellar mass 
densities at high redshift. 

The evolution of stellar mass density will be related to 
that of the star formation rate density in the universe. 
Star formation rate indicator s at the UV wavele ngths are 
sensitive to dust extinction. ICole et all |)2001D adopted 
two cases; E{B - V) = and E(B - V) = 0.15. For each 
case, they provided an analytic formula fitted to the UV 
observations of the star formation rate density. By inte- 
grating these formulae with time, we estimated the evo- 
lution of stellar mass density in the universe. Here, the 
recycling fraction of stellar mass for the next generation 
stars was assumed to be R = 0.28 for a Salpeter IMF. 
The results are overlaid with dashed lines in Fig lA13l 
Lower and upper lines correspond to the E(B — V) = 
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and E(B — V) = 0.15 cases, respectively. It should be 
noted that these lines are shown with units of M Q Mpc -3 
with no dependence on the Hubble parameter because of 
the cancellation in the time integration. 

These time-integrated values of the star formation rate 
densities are derived by the integration of the full range of 
star formation rate or luminosity at each epoch. When 
we assume E(B — V) = 0.15, the line calculated from 
the time integration of the star formation rate densities 
becomes always higher than any of the stellar mass den- 
sity points estimated for the full or the quasi-full range 
of stellar mass. This might suggest that the mean dust 
extinction value would be lower than E[B — V) = 0.15. 
In fact, a median E(B — V) value for the z ~ 4 Lyman 
break galaxies is so mewhat lower than E(B — V) = 0.15 
ijSteidel et alJ ll999') . However, it should be noted that 
all the full or the quasi-full integrated values for the stel- 
lar mass densities at high redshifts were derived from 
the rest-frame UV or optical light only. In addition 
to underestimation of this light due to a certain level 
of dust extinction, there remains some possibility that 
highly reddened systems were completely neglected due 
to their non-detections. In fact, a star formation rate 
density estimate comparable to the values for the case of 
E(B — V) = 0.1 5 was obtained with only a few submil- 
limeter sources (jHuehes et al.lll9 98f). No contributions 
from such dusty galaxies were taken into account in the 
iCole et aTl (|2001^l formulae. It should be noted that our 
stellar mass density estimates include contributions from 
the three submillimeter galaxies in this field. At 6.7/im, 
we can probe the rest-frame near-infrared light out to 
quite a high redshift. Because dust extinction at the 
rest-frame near-infrared is almost negligible, underesti- 
mates of stellar mass densities are unlikely to happen 
either by loss of light or by non-detections. 

6. DISCUSSION 

The contributions of the faint 6.7 fim galaxies to the 
stellar mass density in the universe are estimated to 
be comparable to those inferred from observations of 
UV bright galaxies. Unfortunately, 6.7/im observations 
probe a narrow mass range and UV observations suffer 
from dust effects, preventing us from making detailed 
comparisons between the two. 

On the other hand, we found that the faint 6.7/im 
galaxies generally had red colors. A comparison with a 
particular population synthesis model suggests that they 
have experienced vigorous star formation at high red- 
shifts. The derived large stellar masses for the faint 
6.7 /jm galaxies also support such star forming events 
at the past. Beyond the redshift range of our sample 
[z > 3), we know that there exist Lyman break galaxies. 
They are blue and forming stars; however, their masses 
are generally smaller than the masses of this faint 6.7 /jm 
sample. In a naive sense, several Lyman break galaxies 
must merge to form a massive 6.7 /xm galaxy. Or, very 
rapid star forming systems are needed. SCUBA galaxies 
are expected to have such efficient star forming activi- 
ties. However, at least in this field, SCUBA sources were 
already identified as faint 6.7/im galaxies at relatively 
small redshifts. 

We noticed the existence of massive galaxies out to z = 
3. But at the same time, their comoving space densities 
were found to be lower than the present values. Thus, 



there should be some mass assembly activities to build 
up massive galaxies in a redshift range of z = 0-3. They 
will be forming stars in situ, merging, or a combination 
of the two. The investigation of such build-ups in z — 
0-3 would give us an important insights to the build- 
ups in the higher redshift regime. A large number of 
star forming galaxies detected at 15 /im with ISO may 
give an unbiased sample for this purpose because of their 
insensitivity to dust. 

The detection of all the X-ray, submillimeter, and ra- 
dio sources in this field at 6.7 /im is interesting. Most 
of our knowledge of the evolution of galaxies has been 
based on investigation of stellar systems, such as opti- 
cal observations of UV emission from massive stars. But 
a higher fraction of active galaxies detectable at other 
wavelengths in the distant universe requires the con- 
sistent understanding of the evolution of multiple com- 
ponents in galaxies. Far-infrared/submillimeter obser- 
vations of dust and X-ray /radio observations of active 
galactic nuclei (AGN) might provide us with more essen- 
tial information than UV/optical/near-infrared observa- 
tions for stellar components. 

We have stated that mid-infrared observations re- 
sulted in an investigation of a narrow range of stel- 
lar masses. We regard this as a good property, i.e., 
allowing the preparation of a mass-ordered sample to 
compare stars/dust/ AGN with multi- wavelength obser- 
vations. With mid-infrared surveys in the near future, 
we expect to construct well-controlled samples of distant 
galaxies. 

7. CONCLUSIONS 

The tight correlation between stellar masses of galaxies 
and their rest-frame near-infrared luminosities can be a 
strong tool for investigating the evolution of stellar mass 
assembly in galaxies. The effect of redshift has motivated 
us to observe high redshift galaxies in the mid-infrared. 
For the mid-infrared sources detected in the SSA13 field 
with the ISOCAM LW2 (6.7/im) filter, their nature and 
stellar masses were discussed after their identifications. 
The 65 sources were divided into two subsamples, a pri- 
mary sample of 33 sources (of which 2 are stars) and a 
supplementary sample of 32 sources. No spurious sources 
are expected in the primary sample. Taking account of 
higher source densities at the optical and near-infrared, 
the identifications at these wavelengths were determined 
by utilizing two possibilities; one for the true association 
and the other for a chance event. Using the highest ratio 
of the two, 32 out of the 33 primary sources and 25 out 
of the 32 supplementary sources were identified. A test 
of the identification procedure with the negative sample 
showed that for the supplementary sample identifications 
at the K band should be secure, while a few identifica- 
tions at the / band might be wrong. A comparison with 
the published X-ray, submillimeter, and radio catalogues 
resulted in mid-infrared identifications of all (four) X- 
ray, (three) submillimeter, and (one) radio sources in the 
field. They were all in the primary sample. 

With the color information from the optical and near- 
infrared identifications, we can divide 6.7/im galaxies 
into three types. The GRASIL galaxy model predicts 
that red B — I colors and blue K — 6.7/im colors can be 
used to isolate low redshift early type galaxies (type I). 
Red B — K colors and red I — 6.7/im colors were used 
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to select high redshift early type progenitors (type II), 
which can be regarded as ancestors of the type Is. Blue 
B — K colors and blue 7 — 6.7 /im colors are an indicator 
of on-going star formation (type III). The main contribu- 
tors to the type III category would be late type galaxies. 
The ratio of the three type was almost 1 : 1 : 1 in the 31 
primary galaxies, while the supplementary galaxies have 
a higher fraction of type II. 

In order to permit quantitative discussions, we es- 
timated photometric redshifts for the 6.7 /jm galaxies. 
Based on their good representations in the color-color 
plots, we used a limited set of SEDs in the GRASIL li- 
brary for the fitting templates. A test of this photometric 
redshift estimation with a spectroscopic subsample sug- 
gested that they could be used as good redshift estimates. 
Although a few of them could be high redshift outliers, 
the large errors in their photometric redshifts might be 
used as indicators of such. 

Modulo the caveats on photometric redshifts, we found 
that deep 6.7 /im surveys were efficient in detecting high 
redshift galaxies. The flux redshift relation of the pri- 
mary 6.7 /im galaxies showed a high redshift tail at fluxes 
below 30/iJy. A significant fraction of the photomet- 
ric redshifts for the supplementary galaxies were z > 1, 
consistent with the dominance of type II in them. A 
K < 20 sample has low redshift galaxies with K > 18 
and I — K < 2, while our 6.7 /xm sample do not include 
such a population of blue low-mass galaxies. 

Stellar masses were derived utilizing a tight correlation 
between rest-frame near-infrared luminosity and stellar 
mass. Stellar mass-to-light ratios were determined from 
fits to a limited set of template SEDs in the GRASIL 
library. With a hybrid conversion to stellar mass, using 
K luminosities at z < 1 and 6.7 /im luminosities at z > 1, 
we estimated stellar masses for the 6.7 /im galaxies. We 
found that some of the high redshift 6.7 /im galaxies had 
stellar masses comparable to the typical stellar mass of 
local galaxies. However, the comoving space density of 



such massive galaxies is likely to be a decreasing function 
of redshift. The epoch-dependent stellar mass functions 
might suggest that more massive galaxies are rarer at 
higher redshifts. If our photometric redshifts are slight 
overestimates, this trend will be even stronger. 

We derived the contributions of the 6.7 /im galaxies 
to the stellar mass density in the universe as a func- 
tion of redshift. Given the narrow mass ranges, our es- 
timates were obtained as simple summations of the de- 
tected sources. Our low redshift value was almost con- 
sistent with the local value, suggesting that our sam- 
ple includes major contributors to the low redshift stel- 
lar mass density. At the same time, most of the mass 
assembly in galaxies should be finalized at this epoch 
(z ~ 0.3). The stellar mass density estimates become 
smaller at higher redshifts. This would be expected from 
the decrease in the high redshift stellar mass functions; 
however, the density we compute for the highest red- 
shift bin is almost comparable to the full mass range 
value for a rest-frame optical selected sample. Note that 
our value includes contributions of dusty submillimeter 
galaxies. A full mass range sample of rest-frame near- 
infrared selected galaxies will be necessary to estimate 
correct values of the stellar mass density taking into ac- 
count any dusty population. 
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APPENDIX 

THE GRASIL SED LIBRARY 

Granato and Silva introduced effects of Graphite and Silicate dust to a SED model of galaxies based on 
the stellar population synthesis method <)Silva et alJ 119981) . Effects of dusty interstellar media are treated with 
a radiative transfer code. Thus, the GRASIL SEDs cannot be arbitrarily rescaled. The modeled SEDs in- 
clude UIB emissions, which are believed to be caused by policyclic aromatic hydrocarbon (PAH) particles. A 
small suite of SEDs and executables to calculate SEDs are publicly available at the official GRASIL web site 
(http://web.pd.astro.it/granato/grasil/grasil.html). The full details of the model can be traced at that 
site and the SED library itself is expanding. In order to avoid generating inadequate SEDs by using executables with 
excessive parameter sets, we just have used the author-proofed SED library in this paper. Here we describe the salient 
properties of two sets of SEDs in the library. 

One is a set of evolving SEDs for four Hubble types in the local universe; E, Sa, Sb, and Sc. Elliptical galax- 
ies are assumed to be formed in a mon olithic collapse scenario. The classification of spiral galaxies is based on 
iSolanes. Salvador-Sole, k, Sanromal l)1989|) . These SEDs span from UV to radio wavelengths. The number of ages are 
12 for the E galaxy (0.1, 0.2, 0.4, 0.8, 1.5, 2, 3, 4, 5, 8, 11, 13Gyr) and 15 for the Sa, Sb, and Sc galaxies (1, 2, 3 4, 5 
6, 7, 8, 9, 10, 11, 12, 13, 14, 15Gyr). Their star formation histories are assumed to be smooth, as shown in Fig. IA141 
By setting a cosmology and a formation r edshift, we can convert their age into redshift. 

For stellar populations, Salpeter IMFs i|Salneterlll955[) are assumed with mass limits of 0.15 M Q and 120 M© for 
the E ga laxy and 0. 1 M© and 100 M© for the Sa, Sb, and Sc galaxies. Using a K band luminosity function for local 
galaxies. iCole et"al"l l|2001f) have provided a stellar mass function for a Salpeter IMF with mass limits of 0.1 M© and 
125 Mq. The characteristic stellar mass for this stellar mass function becomes M* tar = 1.7 x 10 11 M© for our adopted 
cosmology. With this unit, the GRASIL E, Sa, Sb, and Sc galaxies with a formation redshift of Zf = 2 have stellar 
masses of ~ 3M* tar , ~ 0.8M* tar , ~ 0.7M* tar , and ~ 0.2M* tax at z = 0, respectively. All the stellar mass estimates 
with the GRASIL SEDs in this paper were normalized to a Salpeter IMF with mass limits of 0.1 M© and 125 M©. 
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The other set of SEDs contains fitted templates of nearby galaxies, which are detailed in iSilva et~aT] llT99l . The 
model parameters were chosen to represent SEDs for starburst galaxies (M82, NGC6090, and Arp220), normal galaxies 
(M51, M100, and NGC6946), and a giant elliptical. The wavelength range is shorter than the evolving SEDs above, 
from UV to submillimeter. 

We used all the SEDs above to derive photometric redshifts. To obtain stellar masses, we only utilized the evolving 
SEDs, since some of the fitted templates lack stellar mass information. Model predictions based on the evolving SEDs 
are shown in some plots, though the Sb SED was omitted to avoid overcrowding. 
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Fig. Al. — Thumbnail im ages of the identified 6.7 /an sources. The images were taken from the data where the actual identifications 
were performed (IDs in Tabl aAlt . This imaging band and the name of the 6.7/^m source were overlaid on each 18 X 18arcsec panel, an 
area of 3 X 3 pixels for the ISOCAM 6 arcsec pixel-field-of-view mode. The identified sources were centered in the images. North is up and 
east to the left. 
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Fig. A2. — f 1/(6-7 /im) / f v (K) vs B — I colors for the identified 6.7 (im galaxies. The sources in the primary and supplementary 
samples (Sect.|5J are shown with filled and empty circles, respectively. One sigma error bars are marked only for the primary sources. For 
sources which have enough information, their color limits are also indicated. We overl aid m odel predictions with the GRASIL SED library 
(Appendix[XJ. Galaxies with star formation histories of the E, Sa, and Sc types CFig. lAlU are shown with dashed, dotted, and solid lines, 
respectively. The assumed formation redshifts are Zf = 2, 3, and 10 (thick, medium, and thin lines, respectively). Triangles mark positions 
for z = 0, 1, 2, and 3 (for z < Zf ; filled triangles for z = 0). In the upper-left region of the plot separated by the dot-dashed line, i.e., 
B — I > 3 log(/„ (6.7 fiTa)/f v (K)/2) + 5, we can identify a group of galaxies marked with squares (type I). They are consistent with the 
evolving E galaxies at z = 0—1, or the Sa galaxies at z = 0-1 with high formation redshifts. It should be noted that the ages of galaxies 
with Zf = 2, 3, and 10 at z = 1 are 2.5 Gyr, 3.6 Gyr, and 5.3 Gyr, respectively, in our adopted cosmology. Thus, type I galaxies will be 
matured systems with old stellar populations created by vigorous star formation episodes a long time ago (see Fig. lA14l . 
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Fig. A3. — fv(6.7 [im)/ fv(I) vs B — K colors for the identified 6.7/im galaxies. Symbols are the same as in Fig. IA2 J alt hough plotted 
sources are slightly different due to different color sets used. This diagram divides galaxies in the lower-right region of Fig. lA2l into two types. 
All the sources above the dot-dashed line, i.e., B — K > —1.7 log(/„(6.7 /im)//» (/)) + 7.3, which have red colors both in /„(6.7/im)//i; (I) 
and B — K, are marked with diamonds (type II). The GRASIL SED predictions show that type II galaxies are z > 1 galaxies with intense 
star formation (evolving E or Sa galaxies). They are consistent with being ancestors of type I galaxies (squares). The remaining blue 
galaxies are categorized as type III. They have the colors of the evolving Sc galaxies or the E or Sa progenitors at 2 ~ Zf, suggesting 
on-going star formation. 
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Fig. A4. — Some examples of the SED fits used to derive photometric redshifts. Detections are shown with squares, in which flux 
errors and band widths arc indicated with vertical and horizontal bars. Arrows are 3 <r upper limits. Source names and their photometric 
redshifts z p h are given with their color types, GRASIL SED names, and \ 2 values (in parentheses). 
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Fig. A5. — Co mpari son of photometric redshift estimates to spectroscopic redshift measurements. The photometric redshifts were 
derived as in Sect. 14.2.11 and their nominal 90 % confidence limits are shown with vertical bars. Almost all the estimates are within one 
sigma ranges (0.2 dex; dotted lines) from the unity relation (solid line). The primary and supplementary sample galaxies (Sect. |5J are 
shown with filled and empty symbols, respectively. 
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Fig. A6. — 6.7 fim fluxes as a function of redshift. Sources with spectroscopic redshifts are indicated with double circles. For photometric 
redshifts, their 90 % confidence limits are indicated with horizontal bars. Vertical bars are one sigma flux errors. Galaxies in the primary 
sample are shown with filled circles, while the supplementary sources are m arked with empty circles with no error bars. Letters, 'X', 'S', 
and 'R' indicate the X-ray, submillimeter, and radio counterparts in Table lA*2l respectively. Dashed, dotted, and solid lines are model 
predictions for the GRASIL evolving E, Sa, and Sc galaxies, respectively (Appendix^}. For each SED, formation redshifts are assumed to 
be Zf = 2, 3, and 10 (thick, medium, and thin lines, respectively). 
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Fig. A7. — K magnitudes of the 6.7 fim galaxies. Symbols are the same as in Fig. IA6l Som e of the 6.7 /an sources are shown with their 
upper o r lower magnitude limits, or not shown due to the lack of their K photometry (Table lA~3l . For reference, a K > 20 spectroscopic 
sample (Cowie et al. 1996) is overlaid with small squares. The 6.7 /xm-selected sample does not include faint K sources at low redshifts, 
some of which are even fainter than the GRASIL evolving Sc galaxy predictions. 
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FlG. AS. — I — K colors of the 6.7/im galaxies. Symbols are the same as in Fig. lA7l The 6.7 /^m-selected sample does not include blue 
I — K sources at low redshifts, some of which are even bluer than the GRASIL evolving Sc galaxy predictions. 
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Fig. A9. — Stellar mass-to-light ratios as a function of redshift. They are estimated for observations at the B, I, K, and ISOCAM LW2 
(6.7 fim) bands. The in-band luminosity at each observing band corresponds to different wavelength with increasing redshift. Thus, for the 
calculation of the ratios, we adopt luminosities in units of the bolometric luminosity of the Sun (lLg o1 = 3.85 X 10 26 [W]). The GRASIL 
models for the evolving E, Sa, Sb, and Sc galaxies are shown with dashed, dotted, dot-dashed, and solid lines, respectively. Two formation 
redshifts are assumed; Z{ = 5 and 10 (thick and thin lines, respectively). 
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FlG. A10. — Stellar masses of the 6.7 fim galaxies. Symbols are the same as in Fig. 

ED The stellar masses are derived from 6.7 fim 
luminosities at z > 1 or K band luminosities otherwise. If no K band luminosities are available, / band luminosities are used instead. Note 
that errors in stellar mass come mainly from errors in photometry, not from errors in stellar mass-to-li ght ratios derived from the SED fit. 
A horizontal dot-dashed line indicates a characteristic stellar mass of local galaxies ICole et alJI2001f> . A dashed line shows stellar mass 
predictions assuming the GRASIL evolving E galaxy with a formation redshift of Zf = 10. They are scaled to the minimum 6.7/im flux 
of the primary sample (12 /^Jy). This curve itself is useful to indicate the redshift dependence of the detection limit especially at the high 
redshift end where the scaling from the GRASIL model is almost unity. However, the scaling factor becomes very small, ~ 0.03 (0.1 M* tar ) 
at z = 0.2 (AppendixF^. So, the uncertainties of the dashed line at the low redshift will be large. Moreover, 6.7/xm fluxes start to be 
affected by dust emission at such low redshifts. Note that stellar masses were derived from K fluxes at Z < 1. 
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Fig. All. — Stellar mass function s at four epochs. For a local stellar mass function, the stepwise maximum likelihood estimates for 
a Salpeter IMF are reproduced from lCole et al] 120011) . The accompanying errors are plotted, though almost all are within symbols for 
values (diamonds). This local stellar mass function is reproduced in three higher redshift panels (dashed lines). In each of these panels, 
our stellar mass function estimates are shown with solid and double circles for the combined (photometric and spectroscopic) and for the 
spectroscopic samples, respectively. Double circles are shifted slightly. Vertical errors are uncertainties in V max and Poisson statistics. 
Horizontal errors are mean fractional errors in stellar mass. 




Fig. A12. — Redshift dependenc e of integrated stel l ar ma ss functions. Lower stellar mass thresholds are those in lDrorv et alj t2001l) . 
Our estimates are shown as in Fig. lAlll Drorv et al. (200l)'s values are represented as square s with vertical err or bars indicating their 
field-to- field fluctuations. Diamonds show local values calculated from a stellar mass function bv lCole et aH J200ltl . 
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Fig. A13. — Stellar mass density in the universe as a function of redshift. The right axis shows densities normalized to the critical 
density of the universe. The contributions of the 6.7 /an galaxies are shown with solid and double circles for the combined and spectroscopic 
samples, respectively. The double circles are plotted at slightly lower redshifts. The horizontal bars represent the redshift ranges of the 
bins and the vertical bars show one sigma errors, taking account of Poisson noise and uncertaintie s in stellar mass and V ma x. Several 
other estimates are overlaid; solid squares are a K band se lected sample with p hotometric redshifts (Fontana ct al. 2003), empty squares 
are an H band se lected sample with photometric redshifts (Dickinson ct al. 2003), t he diamond is a J b and selected spectroscopic sample 
ICole e t al. 2001), em pty circles ar e an / band selected spectroscopic sample (Brinchmann & Ellis 2000), X marks are an R band selected 
spect roscopic sample I Cohen 20(3), an d triangles are an R' band selected photometric sample around a z = 4.7 quasar IGiallongo et al] 
1998). The empty/solid squares and the diamond are obtained from full integration of a Schechter fit to their respective luminosity or 
stellar mass function at each redshift bin. The X marks and the empty circles are quasi-fully integrated values with a finite integration 
range from 10 L* to 1/20 L* , and 10.5 < log(Af s tar[fefi fj 2 Mq]) < 11.6, respectively. The triangles are simply summed values of the detected 
sources like ours. We adopted her mean redshifts for Cohen (2002) and all four cases discussed in Dickinson et al.l 120031) for their error 
bars. Miller-Scalo IMFs were used in IGiallongo et alJ 1199a) . All points plotted here were normalized to our cosmology and a Salpctcr 
IMF with mass limits of 0.1 Mq and 125 Mq. The two dashed curves are deduced by integrating th e star formation rate density in the 
universe, which is derived from the UV luminosity density as a function of redshift I Cole et al. 2001). The upper curve is an extinction 
corrected case for E(B — V) = 0.15, and the lower one has no dust correction. 
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Fig. A14. — Star formation rates (SFRs) for the evolving E, Sa, Sb, and Sc galaxies in the GRASIL library (Appendix [XJ • They are 
shown as a function of age with dashed, dotted, dot-dashed, and solid lines, respectively. It should be noted that age increases from left to 
right. After the age of 1 Gyr, the E galaxy was assumed to stop its star formation because of a galactic wind. 
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Table Al. Source Identifications in the Optical and Near-Infrared — Contin- 
ued 
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Note. — Source identification results with the probability ratio method described in Scct. fTTl The contents are divided into three sub samples; 
the p rimary, sup plementary, and negative samples (Sect.[2j. Names, coordinates, their one sigma errors (ro), and total 6.7/i.m fluxes are extracted 
from ISato et al ] i2003h . The coordinates of the successfully identified sources arc updated with their optical or near-infrared positions in the 
identifying image, which is indicated in the ID column. Listed signal-to-noisc ratios (S/Ns) are derived from standard deviations of the co-added 
images at the source position using a 7.2 arcscc aperture. Thus both photon noise and undersampling effects arc included. With raster observations, 
pixel values could range from ~ 25 % to ~ 100 % of source fluxes depending on exact locations of the sources on the ISOCAM 6 arcsec pixels, giving 
high standard deviations at source positions. At 6.7/xm, an Airy disk size is 5.6 arcsec for the ISO telescope. Note that errors in total 6.7 f-Lm fluxes 
ar e much larger tha n that can be inferred from these detection S/N values. Derivation of total fluxes and their associated errors are described 
in Sato et al For each identification, we list the distance between the 6.7,um source and the identified source in the ID image (ri). Two 

probabilities that the identified source is a true association (Pt) or a chance association (P c ), and their ratio are also shown. The last column shows 
the second highest probability ratio for sources within a distance of 3 r$ from the original 6.7 /im position. The contents for the negative sample 
are the same but with negative flux and S/N values. 



Table A2. Source Identifications at the X-ray, Submillimeter, and Radio Wave- 
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28 


27 










7.7 


J131225. 7+424350 


< 25 


2.4 








40 


14 


0.5 


15 


5.7 


< 4.2 


4.7 


J131228. 0+424458 


< 25 


2.3 








57 


18 










1.8 


,1131232.1+424430 


< 25 


3.8 








62 


88 


0.3 


10 


7.8 


< 3.9 

















Note. — Source identifications of the mid-in frared (ISOCAM/LW2) sources in the SSA13 field with X-ray ( Chand ra) , submillimeter 
(JCMT/SCUBA), and radio (VLA/FIRST) sources IBarger et"allll99a 119991 l200ll IMushotzkv et al.ll2000t IWhite et al.lll997D . No data are in- 
dicated with ' ■ ■ ■ '. Flux limits are n o de tections at a level of 3 <j. Listed offsets dx, <^S; and tZfi are measured from the optical or near-infrared 
positions of the 6.7 /im sources in Table IaTI 
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Name 6.7/^m K I B Type z Bp , z ph SEDi X \ log(M star ) SED 2 xl 

# (/^Jy) (mag) (mag) (mag) (/lg 5 2 M©) 



Primary Sample 



3 50 15.8 18.7 22.1 I 0.314 •■• •■• ll.16jlo.o6 E8 4 - 7 

6 66 16.1 19.0 22.1 I 0.313 •■• •■• 11.08j°°g Sa 10 2.2 

9 25 18.7 22.6 24.7 II 1.31 t° \l Sb 8 0.0 (d4/4) 11.39j°;^ Sa 5 0.6 

10 21 18.8 22.7 > 25.5 II 0.914 •■• •■• 10.89j°°g Sa 6 8.5 

11 19 17.4 20.3 24.6 I 0.732 •■• •■• 11.09j°J| E3 2.4 

13 24 17.3 20.6 24.5 I 0.479 • ■ • • ■ • 10.991° ° i E 8 13.7 

14 19 18.9 22.0 23.6 III 1.305 •■• •■• 11.05j°;^ Sa 2 0.1 

15 60 17.4 20.1 22.3 III 0.279 •■• •■• 10.27j° ™ Sb 6 2.0 

16 19 19.6 22.3 23.9 III 1.168 •■• • ■• 10.96j° ; 22 Sb 3 2.9 

17 17 19.9 24.5 > 25.5 II 2.12jJ;|* Sb 8 0.0 (d3/4) 11.22j°J| E 1.5 1.2 
20 20 19.0 23.0 > 25.5 II 1.16j°;^ Arp 220 0.3 (d3/4) 11.28j°;^ E4 1.8 
22 135 13.8 < 17.0 18.6 star 0.000 

25 39 16.4 19.1 22.4 I 0.24j°°| Sa 15 0.1 (d4/4) 10.75j°°g Sa 11 4.3 

27 45 ■•■ > 23.7 > 25.5 II 

28 27 18.5 21.6 24.3 III 1.038 •■• •■• ll.26jo.2i Sb 6 11.4 
30 70 •■• 18.2 21.4 I 0.277 •■• •■• ll.llj°^ E8 0.1 
33 17 19.2 23.7 > 25.5 II 1.43j{J'.22 Sb 13 0.0 (d3/4) 11.28lo" 2 o E 4 °- 4 
35 15 18.8 22.1 24.6 III 1.08j°°i Sa 4 0.0 (d4/4) 10.98jJ5"i9 Sa 4 
39 27 •■• 22.9 25.4 II 2.96lg 2 | Sa 2 0.0 (d3/3) 11.67j 3o Sa 2 00 

14 19.9 24.0 > 25.5 II 2.70j ®| Sb 4 0.2 (d4/6) 11.32 j°; 2 t Sa 2 4.4 i 

46 16.2 18.9 22.6 I 0.489 •■• •■• 11.33±gJ| E5 2.3 

37 < 17.8 20.0 22.2 III 0.38j° Sc 6 0.0 (d3/4) 10.29j°J Sc 6 0.0 

47 16 17.5 20.5 24.1 I 0.572 •■• •■• 10.9lj° j; E4 2.9 

49 15 19.1 22.8 > 25.5 II 1.05j°;^ Arp 220 0.0 (d3/4) 11.08j°;^ Sa 6 1.0 i 

50 19 18.5 21.4 24.3 III 0.612 •■• •■• 10.70jo°3 Sb 8 0.5 

51 22 17.2 20.4 23.9 I 0.378 • ■ • • ■ • 10.80j° °g Sa 10 9.5 

54 167 15.6 18.4 20.9 III 0.56j°; 2 ^ NGC 6946 0.2 (d4/4) 11.58j° ™ Sa 4 2.7 i 

55 12 18.3 20.9 22.9 III 0.839 • ■ • • ■ • 10.7lj° ™ Sc 4 4.2 
57 18 19.7 > 24.7 > 25.5 II 1.94j°;^ Arp 220 0.0 (d3/6) 11.03 j°;^ E 0.1 6.2 
59 18 17.2 20.1 23.6 I 0.415 •■• •■• 10.88j°°g E8 0.7 

61 55 16.8 19.4 21.5 III 0.276 •■• •■• 10.44j°™ Sa 4 6.6 

62 88 •■• 19.3 22.4 I 0.241 •■• •■• 10.64j°;°* Sail 2.1 
64 107 < 13.7 < 17.2 19.4 star 0.000 



40 
41 
44 



30 



Y. Sato et al. 

Table A3. Properties of the Mid-Infrared Sources — 



Continued 



Name 6.7/^m K I B Typo z Bp , z ph SEDi x \ log(M star ) SED 2 xl 

# OJy) (mag) (mag) (mag) ( h 65 M @) 



Supplementary Sample 

1 
2 
4 
5 
7 
8 
12 
18 
19 
21 
23 
24 
26 
29 
31 
32 
34 
36 
37 
38 
42 
43 
45 
46 
48 
52 
53 
56 
58 
60 
63 



19 
18 
14 
18 
34 
66 
12 
63 
8 
17 



12 
10 
8 
6 
9 
45 
8 
87 
10 
32 
12 
7 
11 



15 
26 
18 



16.7 

> 19.9 

> 19.9 
19.7 



17.9 

19.0 
18.1 

> 19.9 
19.5 

> 19.9 
19.2 

> 19.9 

> 19.9 
19.8 

19.4 

17.6 



18.7 



9 > 19.9 
24 
27 



19.0 



18.4 



19.5 

> 24.7 
24.4 
24.1 

> 23.7 

> 23.7 
21.7 
22.9 
23.6 
21.7 
24.7 
23.4 
25.0 
23.0 
22.6 
22.7 

> 23.7 
23.1 
24.8 

> 23.7 
21.6 

> 23.7 
22.7 
22.5 
23.4 

> 24.7 
23.5 

> 23.7 
22.5 
23.7 
22.2 
23.4 



22.7 

> 25.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 
24.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 
24.2 
25.6 

> 25.5 
24.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 
25.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 

> 25.5 
25.2 

> 25.5 
24.0 

> 25.5 



0.314 



10.84to.o9 Sa 10 3.6 (d4/4) 



2.237 



2.18 



+ 1.92 
0.52 



1.028 



1 93+ 035 



1.06 



+ 0.15 
0.09 



76+ 08 

u. I u_ 08 



1.73 



+ 0.47 
0.52 



1.08! 



9 7n+ 3 - 21 

A - ,u -1.03 



1.02! 



0.620 



0.76^ 



^ + 0.86 
-0.67 



2.09^ 



1.39! 



0.75! 



1 07+ 014 

1 - u ' -0.23 



0.81! 



, 00 + 1.41 
J,M -1.19 



1.31 



+ 0.23 
0.34 



0.88 



+ 0.08 
0.17 



1 43+ - 29 

1-"* J _0.46 



Sa 4 
Sb 6 



I 1.33_°; 2 



, 99 + 0.21 



Sa 7 
GE 
Sa 15 
Sa 7 
Sb 11 
Sb 6 
Sb 11 

Sc 12 
GE 
Sb 6 
GE 

Sa 15 

Sa 7 
GE 
E 0.8 

GE 

M51 
E 5 
Sc 10 
GE 



0.0 (d2/4) 11.27_|| 



0.0 (d3/3) 
0.4 (d3/4) 
0.0 (d3/4) 
0.0 (d2/4) 
0.0 (d3/4) 
0.0 (d2/4) 
0.0 (d3/4) 



+0.18 
27 



0.0 (d3/4) 11.39_ '.2? 



10-93t .20 
H-87ioi 8 o 

10.89to.34 

11 00+ 06 

I l.uu_ 06 

II 01+ 017 

11.U1_Q.28 



10.93! 



ll-18_o'.22 

10 97 +0 18 

lU.y t _ 25 



9.26! 



0.1 (ul/4) 
0.0 (d2/4) 
0.0 (d3/3) 11.76 
0.9 (d3/4) 10.96 



10.28_°™ 



10 97+° 20 

1U.3 I _ 24 



9.9 (d3/4) 11.17_° 



+ 0.05 
09 



Sa 3 
Sa 3 



E 3 
Sa 3 
E 5 
E 5 
E 3 
E 5 
E 1.5 
E 5 
Sb 1 
Sb 6 
Sa 4 
Sa 3 
E 4 

E 5 



Sa 5 

0.9 (d3/4) 10.85_°°5 Sa 7 
0.0 (d2/3) 11.25_°;_^ E 0.8 



0.0 (d2/3) 11.48_°;_| E 5 



0.1 (d4/4) 10.78_ o ° 7 Sa 6 

0.0 (d2/3) II.22I2-30 E 4 

0.0 (d3/3) 11.29+2 22 Sa 3 

0.3 (d3/4) 11.32+°;^ E 5 



0.3 (d2/4) 
1.1 (ul/4) 



3.7 (d3/4) 
2.0 (d3/3) 
3.2 (d3/4) 

4.2 (d3/4) 
0.0 (u2/4) 
0.3 (d3/4) 
0.1 (u2/4) 
0.3 (d3/4) 
0.7 (d3/4) 
0.3 (ul/4) 
0.0 (u2/4) 
0.6 (d3/3) 
4.7 (d3/4) 

16.6 (d3/4) 

0.5 (d3/3) 
4.7 (d3/4) 
0.0 (d2/3) 

0.9 (d2/3) 

1.3 (d4/4) 
0.0 (d2/3) 
0.1 (d3/3) 

5.4 (d3/4) 
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Name 6.7 /im K I B Type z sp , z ph SEDi X \ log(M eta r) SED 2 X 2 

# (MJy) (mag) (mag) (mag) {h~^ M ) 



Note. — Pr operti es of the primary and supplementary 6.7 fim samples. Magnitudes at the K , /, and B bands are corrected total magnitudes as 
in lCowie et al.l 01996ft . Measurements with a 3" aperture are corrected with an average offset. No data arc indicated with ' ■ ■ ■ '. Upper flux limits 
are no detections at a level of 3 a . Lower flux limits are real measurements for saturated sources or for sources at th e im age boundaries. Except two 
stars, the remaining sources are divided into three col or typ es fSect. lTTTl using the GRASIL SED library ( Appendix lAI . Spectroscopic redshifts z sp 
0Barger et al 1 l200lf lC owic et al 1 IT99llSongaii a et al arc listed with no error information, while photometric redshifts z p h are shown with 

their 90% confidence limits. Photometric redshifts z p ^ are derived with fits to the GRASIL evolving templates and local ones [Sect. H:.2.H , while 
stellar masses Af s tar are derived with fits to the GRASIL evolving templates only (Sect. 15.21 . The results of these two fits are shown separately 
as (SEDi, Xi) an d (SED2, X2)- The name of best-fit SED is listed with its age in Gyr for the case of the evolving SED. The local template for 
giant ellipticals is marked as GE. Just after the % 2 value, ancillary information on the fit are shown in parenthesis. For example, d3/4 means that 
3 detections among 4 photometric data were used to determine the best-fit and the resulting x 2 value. The remaining one flux limit was checked 
its consistency to the fit. In case of ul/4, one upper flux limit (non-detection) found to be inconsistent with the initial best-fit derived with the 
remaining three detections. Then, the normalization of the initial fit was decreased down to the upper limit and the x 2 value was evaluated with 
this re-normalizcd spectrum. 



Table A4. Stellar Mass Functions and Stellar Mass 
Densities 



z Range log $(M star , z ) logp sta r(^) 

(/4 Mpc -3 dex" 1 ) (/i 6 5 M Mpc" 3 ) 

logM star range 

10.15 - 10.55 10.55 - 10.95 10.95 - 11.35 11.35 - 11.75 



Combined Sample 

0.2-0.5 -2.3i°; 5 5 (3) -2.2±j>;»(4) -2.3±g;| (5) ■•■ 8.7±g;l (12) 

0.5-1.2 ■ • • -2.5t°; 5 6 (4) -2.9«; 3 3 (6) -3.81?;* (1) 8.2l°j (11) 

1.2-3.0 ■ • • • • • -3.3l<j:5 (5) -4.018;? ( 2 ) 7.6±g;t ( 7) 

Spectroscopic Sample 

0.2-0.5 -2.518:3 (2) -2.418 J (3) -2-3+°;g (5) 8.61° j (10) 

0.5-1.2 ... -2.5l8: 5 6 (4) -3.2181(3) ••• 8.1181(7) 

1.2-3.0 ••• ••• -3.71^(1) ••■ 6.9lS-J(l) 



Note. — Epoch dependent stellar mass functions $(M star , z) and stellar mass densities p s ta,r(z). They are derived 
from the 6.7 /im primary sample. Two sets of estimates are obtained; one with the combined (photometric and 
spectroscopic) sample and the other with the spectroscopic sample. One sigma errors are estimated from uncertainties 
in Vmax and Poisson statistics. For stellar mass densities, errors in stellar mass are also taken into account. Numbers 
in parentheses are the numbers of 6.7 fim galaxies used in the calculation. 



